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2 ' Abstract. In this paper we present a new deep, wide-field near-infrared imaging survey. Our J- and A^-band observations in 

' four separate fields (0226-04, 2217-1-00, 1003-1-02, 1400-H05) complement optical BVRl, ultraviolet and spectroscopic observa- 

, tions undertaken as part of the VIMOS-VUT deep survey (VVDS). In total, our survey spans ~ 400 arcmin^. Our catalogues 

are reliable in all fields to at least K ~ 20.75 and J ~ 21.50 (defined as the magnitude where object contamination is less than 
10% and completeness greater than 90%). 

Taken together these four fields represents a unique combination of depth, wavelength coverage and area. Most importantly, 
5^ ' our survey regions span a broad range of right ascension and declination which allow us to make a robust estimate of the effects 

of cosmic variance. We describe the complete data reduction process from raw observations to the construction of source lists 
and outline a comprehensive series of tests carried out to characterise the reliability of the final catalogues. From simulations 
we determine the completeness function of each final stacked image, and estimate the fraction of spurious sources in each 
magnitude bin. We compare the statistical properties of our catalogues with literature compilations. We find that our J- and 
A'-selected galaxy counts are in good agreement with previously published works, as are our (J-K) versus K colour-magnitude 
diagrams. Stellar number counts extracted from our fields are consistent with a synthetic model of our galaxy. Using the location 
of the stellar locus in colour-magnitude space and the measured field-to-field variation in galaxy number counts we demonstrate 
that the absolute accuracy of our photometric calibration is at the 5% level or better Finally, an investigation of the angular 
clustering of K- selected extended sources in our survey displays the expected scaling behaviour with limiting magnitude, with 
amplitudes in each magnitude bin in broad agreement with literature values. 

In summary, these catalogues will be an excellent tool to investigate the properties of near-infrared selected galaxies, and such 
investigations will be the subject of several articles currently in preparation. 
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1. Introduction 
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2. Observations 



Galaxy surveys selected in near-infrared wavelengths (1.1 - 
1.5/z) provide some well-established advantages with respect 
their optically-selected counterparts. A census conducted at 
these longer wavelengths can provide flux measurements in an 
object's rest-frame optical bandpass at intermediate redshifts, 
which can in turn be easier to relate to physically interest- 
ing quantities such a galaxy's total mass in stars. Any cred- 
ible model of galaxy formation must predict how this stellar 
mass function evolves with redshift. In addition, the predicted 
number densities and spatial distribution of objects lying at the 
reddest outer reaches of the optical-infrared colour-magnitude 
diagram ( "the extremely red objects" or EROs) also depend 
very sensitively on one's assumed model of galaxy formation. 
Understanding this red outlier population and how it relates to 
UV-selected star-forming galaxies at higher redshifts and mas- 
sive ellipticals at the present day has become one of the most 
important questions in observational cosmology. Near-infrared 
data is also crucial to compute accurate photometric redshifts 
in the 1 < z < 2 redshift range, where measuring spectroscopic 
redshifts with a red-optimised spectrograph can be challeng- 
ing. Until recently, however, the small format of near-infrared 
detectors (and the much higher ground-based brightness of 
the sky at longer wavelengths) has made surveys of the near- 
infrared selected Universe a very time-consuming undertaking. 
Only in the last few years, with the advent of larger-format de- 
tectors, has it become practical to survey deeply larger areas of 
the sky in the near IR to cosmologically significant redshifts. 

In this paper we describe a new deep near-infrared sur- 
vey. The observations presented here cover a total area of 
~ 400 arcmin^ over four separate fields in both J and K bands. 
Each of the four fields reaches a completeness limit (defined 
as the magnitude at which 90% of simulated point sources are 
recovered from the images) of at least 22.0 magnitudes in J 
and 20.75 in K. This represents an intermediate regime be- 
tween , for example, very deep surveys like FIRES dLabbe et alJ 
l2003h which covers a few s quare arcminutes to depth of 
K ~ 2 4, and shallower surveys jPaddi et alJbOOOHProrv et alJ 
1200 1 *) which reach to around K ~ 19 over several hundred 
square arcminutes. The current survey has bee n undertaken in 
the co ntext of the VIMOS-VLT deep survey dUe Fevre et alJ 
I2OO5I) . and the near-infrared data presented here comple- 
ments optical and ultr aviolet imagery (McCracken et al. 2003; 
iRadovic hetal] l2004 . as well as 2.4Ghz VLA radio data 
tBondi et al."2003l). To K < 20.5, 766 objects from our cata- 
logue h ave been observed sp ectroscopically in the first epoch 
VVDS dLe Fevre et al.l2005l) . 

Our primary objective in this paper is to describe in detail 
how our near-infrared catalogues were prepared and to quantify 
their reliability and completeness. Future papers will present 
more detailed analysis of our J- and /^-selected samples and 
in particular the clustering properties of objects with extreme 
colours. All magnitudes quoted in this paper are in the Vega 
system unless otherwise stated. 



The survey described in this article covers in J and K 
bands four different regions of the sky spanning a wide 
range of right ascension so that at least one field is ob- 
servable throughout the year. Each of these regions has 
corresponding deep optical imaging, taken with either the 
Canada France Hawaii Telescope's (CFHT) CFH12K camera 
(McCracken et al. 2003; Le Fevre et al. 20 0j) or the E SO 2.2 
metre telescope's Wide-Field Imager dRadovich et alJl2004l) . 
The observations described this paper were carried out at 
the ESO New Techonology Telescope using the SOFI Near 
Infrared imaging camera dMoorwood et alJll998l) with the J 
and the Ks filters. The K, filter is bluer and narrower than the 
standard near-infrared ^T-band filter, and so is less aff'ected by 
the thermal background of the atmosphere and of the telescope 
dWainscoat & Cowiell992l) . Throughout this paper our K mag- 
nitudes are in fact "K/' magnitudes. SOFI is equipped with a 
Rockwell Hawaii HgCdTe 1024x1024 array and observations 
were made with the Large Field (LF) objective, correspond- 
ing to a field-of-view of 4.9 x 4.9 "and a pixels scale of 0.288 
"/pixel. 

The observations took place over a series of runs from 
September 1998 to November 2002. Each targetted field was 
observed in a series of pointings in a raster configuration, each 
separated from suiTounding ones by 4' 15 "in both right ascen- 
sion and in declination, in order to ensure a non negligible over- 
lap between adjacent pointings. Figure shows the layout of 
our A'-band observations of field 0226 - 04 as an example. 

The well known peculiarities of infrared observations 
(higher and more variable sky background, strong and variable 
absorption bands) dictated our observational strategy. Total 
integration time per pointing was around 1 hour for the 
7-band and three hours for /T-band exposures, with some 
pointings being observed for up to four hours in /iT-band. Each 
integration consisted of many shorter jittered exposures with 
the telescope being offset by random amounts within a box 
of size 30". The jittered observations were usually grouped 
in observation sequences each typically one hour long. Each 
individual short exposure in both bands was 1.5 minutes long, 
with DIT = 15 (meaning a detector integration time of 15 
seconds was used) and NDIT = 6 in 7 (indicating that six of 
these 15 second integrations were used) and with DIT = 10 
sec and NDIT = 9 in K. For the observations of the standard 
stars we adopted DIT = 2 ( DIT = 1.2) sec in J (in K) and 
NDIT = 1 5 to avoid saturation. 



The four observed areas are listed in Tabled together with 
the centres for each of the pointings (J2000), observing runs 
when the observations were performed, total exposure times 
and galactic extinction in J - IK- band com puted using the 
COBE/DIRBE dust maps dSchlegel etalJl998l) . Unfortunately, 
poor weather conditions partially hampered our efficiency and 
reduced our final total areal coverage. 
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Table 1. List of fields observed. See text for more details. 



3. Data reduction 

3.1. Science frames 

Data reduction of the scientific frames in both filters included 
the usual standard steps: dark subtraction, flat fielding and sky 
subtraction. We now outhne the basic processing steps fol- 
lowed. The darks to be subtracted were computed using the 
IRAF' task darkcombine from a series of darks obtained with 
the same DIT and NDIT as the science frames. The well known 
complex bias behaviour of the Rockwell Hawaii array means 
that there is a dependence of the detector dark on time and il- 
lumination history, and this manifests itself as a pattern which 
remains in all the images after dark subtraction. This is visible 
as a discontinuity between the lower upper part of the array and 
the upper lower part of the array {i.e. were the two upper quad- 
rants join the two lower quadrants). This pattern is a purely ad- 
ditive component and it is removed in the sky subtraction step 
as it changes little in images in the same stack of observations. 

Flat-field frames were obtained for both bands with the 
ON-OFF procedure as described in the SOFI manual. We care- 
fully checked whether our flat-field frames contained the large 
jump between the lower and upper part of the array which is a 
signature of a remaining dark residual, rejecting those frames 
which had a non-negligible dark residual. Using IRAF's flat- 

' IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Fig. 1. Layout of the observed ^T-band pointings for field 
0226 - 04. The raster configuration chosen in order to secure 
a non-negUgible overlap between adjacent pointings is clearly 
visible. Each pointing is is turn observed through a series of 
jittered exposures, with jitter box size of 30 "(see text for more 
details). 
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combine task we were able to derived from the remaining flats 
the final J or K flat. To test the accuracy of the ON-OFF flat 
field correction, we compared the counts observed for each 
photometric standard star in different array positions. The rms 
on the mean was always of the order of a few percent for both 
bands (and usually below 4%). 

Before proceeding with the data reduction, we checked the 
quality of the images and rejected those with full-width half- 
maximum (FWHM) values larger than l."3 or those which 
were badly affected by sky-transparency fluctuations based on 
the flux measurements of a chosen reference star The exposure 
times quoted in TableQlare those obtained after this step. 

We next used IRAF DIMSUM^ package to subtract the 
sky from the dark subtracted, flat fielded, science images 
llStanford et al.lll995h . DIMSUM was used in a two step pro- 
cess and on sets of 20-30 images belonging to the same obser- 
vation sequence. 

As a first step for each image 6 neighbouring images (3 
minimum for the images at the end/beginning) from the same 
observation sequence were selected to obtain an initial sky es- 
timate. The sky subtracted images were then used to identify 
the image regions covered by the objects. In this way for each 
image an object mask was defined and was used to exclude 
those regions from a second, final pass of sky estimation and 
subsequent subtraction. 

Finally, the second pass sky subtracted images of the 
same observation sequence were combined in a stacked im- 
age by using the task dithercubemean from the IRDR software 
(Sabbey et al. 2001). This task uses bi-linear interpolation to 
register the input frames, taking into account pixel weights, 
based on image variance, exposure time and pixel gain (which 
is calculated using superflat field irnage calculated from all the 
sky subtracted images; see ISabbev et alJ ( l200ll) for more de- 
tails). 

A final weight map was also produced by combining indi- 
vidual image weights. At the end of these data reduction steps 
usually one or two (three or four for K-hwd observations) 
coadded images and their relative weight maps were produced 
for each pointing in /-band, depending on how the observa- 
tions were sequenced. 

Figure 13 shows the FWHM seeing distribution for the dif- 
ferent coadded /T-band images for each of the four fields. 
Tables |2] and |3l list for each field and each band the median 
FWHM seeing as measured in the final mosaicked image (for 
the A'-band images this value is shown by the dashed lines in 
Figure 13 ■ 

3.2. Astrometric calibration 

Astrometric calibration on the coadded images was performed 
in two steps for all our images. We first computed a linear as- 
trometric solution using the astrometric catalogue of the United 
States Naval Observatory (USNO)-A2.0 (Monet 1998) , which 
provides the positions of 0.5 x 10^ sources. The area covered 
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^ DIMSUM is the Deep Infrared Mosaicing Software package, P. 
Eisenhard, M. Dickinson, S. A. Stanford, and J. Ward, available at 
^ftp://iraf.noao.edu/iraf/contrib/dimsumV2/_ 
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Fig. 2. FWHM seeing distribution of the coadded images for 
each field observed in ^T-band. The dashed line shows the me- 
dian seeing measured in the final mosaicked image. Note that in 
all cases this value is below 1". All images with seeing higher 
than l."3 were discarded, and are not shown in this plot. 



by each coadded image usually contained around 10 objects 
(after removal of saturated and extended sources). This first as- 
trometric solution was then improved by using a catalogue of 
sources extracted from the resampled /-band images which is 
used as a referen ce catalogue for all the other optical bands (see 
(iMcCracken et al. , 2003). and the simi lar procedure adopted for 
the f/-band in iRadovich et alJ ll2004l) '). As the surface density 
of these objects is much higher than that of USNO-A2.0, a 
much higher accuracy in the relative astrometric solution can 
be obtained. Such accuracy allows us to match sources at the 
sub-pixel level between optical and infra-red bands. 

The quality of the astrometry in the final /T-band mo- 
saicked images is shown in figure|3]for the four fields. Radial 
residuals between /T-band and /-band positions for unsatu- 
rated, point-like sources are plotted for each of the four mo- 
saicked fields. The inner circle and the outer circle show in 
each field enclose 68% and 90% of all objects, indicating that 
we have reached the level of sub-pixel accuracy in the relative 
astrometry between /-band and /T-band images (the achieved 
RMS positional accuracy, defined as the radius enclosing 68% 
of the objects, is for all four fields around ~ 0."2). Similar 
values of radial residuals are obtained between /-band and 
/-band positions. The RMS accuracy of our absolute astromet- 
ric solution is therefore of the same order of that obtained for 
/-band data, that is ~ 0."3 (see McCracken et al.. (.2003i) V 
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Fig. 3. Radial residuals between /T-band and /-band positions 
for unsaturated, point-like sources for each of the four fields. 
The inner and outer circles enclose 68% and 90% of all objects 
respectively. The dashed lines cross each other at the position 
of the centroid of the residuals distribution. 



3.3. Photometric calibration 

The photometric calibration was performed using standard 
stars from the list of near-infrared NICMOS standard stars 
(Persson et al. 1998). At least two standard stars were observed 
each night at low airmass {sec(z) < 1 .3), at intervals of roughly 
two hours, in both J and K bands. Each standard star was ob- 
served in 5 different array positions: once near the center of 
the array and once in each of the four quadrants. These images 
were dark subtracted and flat fielded according to the procedure 
discussed above. Sky subtraction was obtained by subtracting 
the median of the four adjacent images (usually the standard 
star fields are empty of bright stars). 

Instrumental aperture magnitudes for the standard stars 
were computed within an 8" radius and were airmass cor- 
rected assuming an atmospheric extinction coefficient of 0.1 
(0.05) in magnitudes/airmass in J (K) bands (based on 
recent measurements in LaSilla, see SOFI home page at 
ESO: www.ls.eso.org/lasilla/sciops/ntt/sofi/index.html). We es- 
timated the actual zero point for each night by comparing the 
instrumental magnitudes with those quoted in the literature. 
This way we obtained for each coadded image an absolute pho- 
tometric calibration. 

Not all our nights were of excellent photometric quality, 
and therefore a refinement of this first photometric solution was 
needed. After correcting each coadded image to one second 
exposure time, to zero airmass and scaling to an arbitrary zero 
point (normally we chose 30) we improved on this first solution 
in two steps. 



Fig. 4. Comparison between the magnitude of corresponding 
objects in two different ^T-band coadded images of the same 
pointing. The sequence of images n4 81 - 120 was observed 
during the run of November 1998, while n4 121 - 160 was 
observed during the run of November 2002. The agreement 
for the photometry of these two coadded images is quite good 
(the mean value of the magnitude difference is only Amags - 
-0.015). 

For different coadded images of the same point- 
ing, when available, we compared the magnitude of the 
point-like brighter obje cts in common using SExtractor's 
teertin & Arnoutsl ll996) MAG_AUTO measurements to perform 
the comparison. In this way were able to check for possible 
photometric offsets between different coadded images of the 
same pointing. In case of discrepancies we anchored the zero 
point to the coadded image with the best photometric quality 
(based on the quaUty both of its night of observation and of the 
specific sequence used to build it). Figure 0] shows as an ex- 
ample the comparison between two different coadded /T-band 
images of pointing n4 of field 2217 -f 00, the first taken dur- 
ing the run of November 1998 and the second during the run 
of November 2002. In this case the agreement between these 
two coadded images is quite good. In other cases a shift had 
to be introduced in order to put the different coadded images 
of the same pointing on the same zero point scale. The size of 
these shifts was in a few extreme cases as large as 0.3 mags, 
but mostly below 0.1 mags. 

As a second step, we used the stars in common in the over- 
lapping areas between different pointings to define a common 
photometric scale for each field, assuring the homogeneity of 
the survey photometry. In this case the few brightest non satu- 
rated stars in common enabled us to define a common photo- 
metric solution for the different pointing through scaling fac- 
tors, to be used when building the final mosaicked image of the 
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3.4. Preparation of final staci<s 
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Fig. 5. Root-mean-square magnitude errors as a function of 
A'-band magnitude. See text for more details. 



field. The corresponding shifts introduced in the zero point of 
the fields to be corrected were always below 0.1 mags. 

We took advantage of the repeated observations of the same 
pointing to estimate the random photometric errors in the final 
mosaicked images. Following the same procedure outlined in 
detail in the following Sections, we produced for field 1400-1-05 
3 mosaicked A'-band images, each corresponding to one hour 
of exposure time, and extracted a photometric catalogue from 
each of these images. By comparing the magnitudes measured 
for the same object in the three independent one hour stacks, 
we obtained a direct estimate of the random photometric er- 
rors present in our /T-band data, including flat-fielding and/or 
background subtraction inaccuracies (all the direct error esti- 
mate have been divided by V3 to take into account the shorter 
exposure time for the individual stacks). Figure |5l shows the 
comparison of such direct error estimate with the errors ob- 
tained by SExtractor and with the errors computed from the 
simulations used in Sect. l4.ll to estimate completeness of our 
fields. The errors obtained by SExtractor are always lower, 
usually by a factor of 2, than the more realistic direct error es- 
timate. The errors obtained from simulations are slightly lower 
than the direct error estimate: they refer to stellar objects, and 
are therefore less affected by errors in background determina- 
tion, especially at bright magnitudes. A similar trend is ob- 
served for 7-band magnitudes errors: the estimate provided 
by SExtractor is lower by a factor of roughly two than the 
directly measured errors. A realistic estimate of random photo- 
metric errors as a function of object magnitude has to be taken 
into account when e.g. using these data for photometric redshift 
determination (see Bolzonella et al. 2005, in preparation). 



Once the astrometric and photometric solutions were com- 
puted for all the coadded images of each field, these images, 
together with their weight maps, were combined to produce 
the final stacks and their weight maps. This is carried out in 
a two-step process by using Swarp, an image resampling tool 
(Bertin et al. 2002). This process does not differ from that de- 
sc ribed in detail for the s tacking of optical images described 
in lMcCracken et all ll2003 ). and therefore the interested reader 
should check there for its details. The final mosaicked images 
has the same pixel sc ale (0".205/pixel) and orientation as the 
final optical images jMcCracken et alJl2003 ). For each field 
and band these images were corrected for the mean galactic 
extinction at the cent re of the field as derived from the maps of 
ISchlegel etai1(ll99Rl) . 



3.5. Cataiogue preparation 



We used SExtractor dBertin & ArnoutsI 1 19961) to extract 
sources from the stacked images and their weight maps. For ob- 
jects to be included in our J- or /T-band catalogues they must 
contain at least 3 contiguous pixels above the detection thresh- 
old of 1 . 1 cr, giving a minimum signal to noise ratio per source 
of 1 .9cr for this per-pixel threshold. This conservative threshold 
means we minimise the number of spurious detections while 
not adversely affecting our completeness (see later). For the 
mosaicked image of field 0226 - 04 a ch i-squared BVRIK im- 
age was constructed iSzalavet alJll999h . and for the J- and 
^T-band catalogues, image detection was performed using such 
an image in SExtractor double-image detection mode, with 
similar extraction parameter s as those quoted above (the inter- 
ested reader should refer to ('McCracken et al."2003') for more 
details on the chi-squared image construction). In our cata- 
logues magnitudes were measured using the SExtractor pa- 
rameter riAG_AUTO. This parameter is intended to give a pre- 
cise estimate of total magnitudes for extended objects, an d 
is inspired by Kron's first moment algorithm jKronl 1198(1) . 
We adopted a minimum Kron radius = 1".5, that is 

when objects are faint and unresolved MAG_AUTO magnitudes 
revert to simple aperture magnitudes. Using the same simula- 
tions adopted to test the completeness of our catalogues (see 
Sect. 14.11 ) we verified that MAG_AUTO was a reliable estimate 
of the input total magnitude of the simulated objects and that 
the systematic loss of flux was always smaller, in the range of 
magnitudes of interest, than the dispersion in the magnitudes 
recovered. The catalogues were visually inspected and noisy 
border regions of the mosaicked images were masked out to- 
gether with circular areas surrounding bright objects. 

The total area of our survey, after bad regions are excised, is 
389 arcmin^ in 7-band and 430 arcmin^ in /T-band. The total 
number of detected sources is 6433 down to 7 = 22.00 (7823 
down to 7 = 22.25) and 8105 down to K ^ 20.75 (9539 down 
to K - 21.00). Tables|2]and|3lshow for each field the final area 
covered and the number of objects detected down to different 
magnitude limits. 
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Fig. 6. Half-light radius (ri/2) as a function of /T-band mag- 
nitudes for each of our four fields. Heavy dots show point like 
sources. The histogram on the right hand side shows clearly the 
locus of point-like sources for objects brighter than the classi- 
fication limit. Fainter than this magnitude (K ~ 17) the star- 
galaxy separation becomes unreliable. 

3.6. Star-galaxy separation 

The separation of extended from point-like sources was per- 
formed separately for each field and each band using the 
SExtractor £lux_radius parameter This parameter, de- 
noted as ri/2, measures the radius that encloses 50% of the 
object's total flux. For point-like sources ri/2 is independent 
of magnitude and depends only of the image seeing. As the 
seeing is quite uniform across our mosaicked images for each 
field, a plot of ri/2 vs. J-, or /T-band magnitude clearly defines 
the stellar locus. Heavy dots in figure |6] show, for each field 
in A'-band, the point-like objects selected using this classifier. 
The histogram on the right hand side of each panel illustrates 
for each field the distribution of ri/2 for magnitudes brighter 
than those deemed feasible for a reliable star-galaxy separa- 
tion. The stellar locus is clearly visible in this plot. A similar 
plot was also used to select stars in the /-band images. 



4. Data quality assessment 

In these subsections we present a series of quality assessment 
tests carried out on the catalogues prepared in the previous sec- 
tions. 



4.1. Completeness and contamination 

A simple estimate of our limiting magnitude can be obtained 
by using the background RMS cr provided by SExtractor 
to compute, for each stacked image, the corresponding 3a-- 



(5cr-) magnitude limits. The formula is magincr) - Zp — 
2.5log(ncr yfA), were n=3(5), Zp is the zero point and A is 
the area of an aperture whose radius is the average FWHM 
of point like sources (see Tables |2 and 1^. The values esti- 
mated this way do not vary much from field to field, and are 
mag(3cr) ~ 23. A and magiScr) ~ 22.9 for the 7-band images, 
while are mag{3cr) ~ 22.6 and mag{5cr) ~ 22.1 for the /T-band 
images. These values can be regarded as indicative lower limits 
on the detectability of objects in our catalogues. 

To better characterise the photometric properties of our im- 
ages we carried out an extensive set of simulations. For each 
mosaicked image and for each filter a list of random positions 
was generated and cross-checked with the position of the real 
sources to reject cases of possible overlapping between gener- 
ated random positions and bright real sources positions. The 
remaining list of coordinates was used to add to the stacks arti- 
ficial stars distributed uniformly in the range 17.5 < J/K < 25. 
SExtractor was then run on these images with the same pa- 
rameters adopted for the detection of real objects and the re- 
sulting catalogues were cross-correlated with the input list of 
artificial stars. The process was repeated until a robust statistic 
was obtained (more than 5000 objects per half-magnitude bin) 
and the ratio ninput /"output was estimated in bins on 0.5 mags. 
The results are plotted in Figures and |8l for each field, and 
the values magcomp and magiim, corresponding respectively to 
90% and 50% completeness level, are shown in Tables |2 and 
13 Such completeness limits and curves assume that the pro- 
file of the source is point-like and therefore should be consid- 
ered only as upper limits. The true limiting magnitude of an 
extended source will depend on its light-profile and actual size 
and can be up to one magnitude brighter for low surface bright- 
ness objects, see eg. Cristobal-Hornillos et al.. C2003) . For field 
0226 - 04, where the extraction was done on the BVRIK chi- 
squared image, the values quoted for completeness, being ob- 
tained from the J and /T-band mosaicked images only, are a 
conservative estimate, irrespective of color, for the detectabil- 
ity of objects in the chi-squared image. For each field we also 
measured incompleteness as a function of position across the 
mosaicked image. Such tests did not show any significant vari- 
ation at m = wflgcomp across the mosaicked images (except for 
the small areas around bright stars, those masked out in the fi- 
nal catalogue). 

To complete the characterisation of the photometric prop- 
erties of our fields we estimated contamination rates as a func- 
tion of magnitude bins. Each mosaicked image was multi- 
plied by -1 and detection of (obviously spurious) objects was 
performed on this inverse image using the same parameters 
adopted for detecting objects on the mosaicked image. For 
the field F02, where the detection was carried out using the 
chi-squared image, we used the same technique described in 
iMcCracken et al. (.2003..) . that is the detection of fake objects 
was performed on the inverse image using the chi-squared 
BVRIK image as the reference image. Tables |2 and |3 list, 
for each field and band, mcont<i{)%: the center of the faintest 
0.5 mags bin where contamination by spurious sources is below 
10%. For the magnitude limit for the scientific analysis of each 
of our fields we adopted the conservative choice of the bright- 
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Fig. 7. Ratio of detected input sources as a function of input Fig. 8. Ratio of detected input sources as a function of input 

magnitude for /-band. The input catalogue consists of a flat magnitude for K-hwd. The input catalogue consists of a flat 

distribution of simulated point-like sources. The dotted lines distribution of simulated point-like sources. The dotted lines 

show the adopted 50% and 90% completeness magnitudes for show the adopted 50% and 90% completeness magnitudes for 



each field. 



each field. 



est between mcomp and mcont<i{)%- As expected, field 0226 - 04 
has the lowest contamination rate at the completeness limits 
of the survey, confirming the effectiveness of the chi-squared 
detection image technique in reducing the number of spurious 
detections. 

4.2. Galaxy and star number counts 

Comparing number counts of galaxies and stars with published 
compilations is a good check both of the star-galaxy separation 
efficiency and of the reliability of our photometry, as well as the 
sample reliability and completeness. The diff'erential number 
counts of stars (number 0.5 mag"' deg"^) for each of our fields 
are shown in Figures|5]and^| To avoid underestimating bright 
stars counts, for this exercise we used the catalogues before 
excising the areas around brig ht stars. The continu ous lines are 
the prediction of the model of lRobin et alJ ll2003l) computed at 
the galactic latitude appropriate for each field. The agreement 
between observed and predicted star counts is very good for 
both J and K bands (the error bars shown are Poissonian error 
bars), confirming the reliability both of our photometry and of 
our star-galaxy separation procedure. It should be noted that 
for field 2217 H- 00 the expected star counts are quite high, due 
to the relatively low galactic latitude of this field. For fields 
1003 -H 01 and 1400-1-05 the star counts are slightly lower, while 
for 0226 - 04 we can safely assume that the star contamination 
is well below a few percent for all the relevant magnitude bins 
(see later). 



Figures [H] and [21 show the diff^erential number counts 
(number (0.5 mags)"' deg"^) for our four fields in J and 
K bands, obtained by normalising the observed raw number 
counts to the areas listed in Tables |2 and |3 The error bars 
shown are Poissonian eiTor bars and no coiTection for stel- 
lar contamination has been applied the counts shown. For the 
0226 - 04 field the cont amination estimated from the predic- 
tion of the model of .Robin et alJ (i2003.) is below 5% for the 
fainter bins shown in the plot, while for the brighter ones (be- 
low J = 19.5) rises beyond 10%. The situation is not so favor- 
able for fields 1400-1-05 and 1003+01: for these two fields only 
at magnitudes fainter than J - 20.5 contamination rates falls 
well below 10%. The worst case is field 2217 + 00 where, due 
to the lower galactic latitude of this field only at magnitudes 
fainter than J - 21.0 do contamination rates become negligi- 
ble. A similar trend holds for the ^T-band counts. Having taken 
into account these eff'ects, the agreement among the fields is 
quite remarkable. The dotted line shows, both in FigurefTTIand 
in Figure fT2l the final raw counts (number 0.5mag"' deg"^), ob- 
tained by simply adding up counts from the four fields, while 
the heavy continuous line shows the total galaxy counts ob- 
tained after correcting each field counts for stellar contam ina- 
tion ac cording to the predictions of the model of Robi n et alJ 
ll2003h. Giv e n the excellent agreement between the model of 
'Robin et alj ( l2003l) and our bright star counts (see Figures |5] 
and not , we always use this model to coiTect for stellar con- 
tamination, even in the brightest magnitude bins. On the right 
of the continuous line, slightly offset for sake of clarity, are 
shown +/ - \ (T error bars for total galaxy counts in each bin, 
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field name seeing mcomp mum fncom<\o% Area Ntot N(m < 22.00) N(m < 22.25) 
arcsec arcmin- 



0226 - 04 


1.05 


22.25 


23.25 


23.00 


162.6 


12327 


2604 


3144 


2217 + 00 


1.00 


21.90 


22.90 


21.75 


102.5 


6676 


1659 


1949 


1003 + 01 


0.90 


22.25 


22.90 


21.50 


102.7 


8072 


1689 


2125 


1400 + 05 


1.20 


22.00 


23.15 


21.90 


23.4 


1589 


336 


408 



Table 2. For each field for /-band data this table lists magcomp^ i-^- the center of the 0.5 mag bin where 90% of point-like sources 
are retrieved, magu^ i.e. the center of the 0.5 mag bin where 50% of sources are retrieved, and 'W(,ont<io%^ i-^- is the center of the 
0.5 mag bin where contamination by spurious sources is below 10%. See text for more details. 



field name 


seeing 
arcsec 




"llim 


mcont<lO% 


Area 
arcmin- 


Ntot 


N(m < 20.75) 


N(m < 21.00) 


0226 - 04 


0.82 


20.75 


22.00 


21.50 


168.3 


12755 


3161 


3757 


2217 + 00 


0.87 


21.00 


22.15 


20.75 


88.3 


5850 


1838 


2113 


1003 + 01 


0.80 


20.75 


21.90 


20.50 


122.0 


9114 


2235 


2654 


1400 + 05 


0.73 


20.75 


21.90 


20.50 


45.0 


3486 


733 


870 



Table 3. As for Table|2but for /T-band data. See text for more details. 
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Fig. 9. Differential number counts of stars in 7-band for our Fig. 10. Differential number counts of stars in /T-band for our 
fields. The continuous lines are the prediction of the model of fields. The continu ous lines are the prediction of the model of 



[Robin et al. ( 2003) for each field. 



iRobin et alJ ( l2003l) for each field. 



obtained by computing the weighted variance among corrected 
galaxy counts for each field (with weighting proportional to 
area covered). The size of these error bars indicates that taking 
into account the different areas covered for each field, and after 
correcting for stellar contamination, results from each field are 
consistent at roughly the 5% level. We did not need to apply 
any other completeness or contamination correction to our data 
because down to the limiting magnitude plotted in Figures [O] 
and[^such corrections are negligible. 



Tables ^ and [S] list our raw differential number counts in 
each field, the total raw number densities (in units of num- 
ber 0.5 mag ' deg"^), and the final, corrected for stellar con- 
tamination, galaxy densities together with their Icr error bars, 
computed as described above, for our total sample. 

Figures [T3l and [T4l show our total corrected galaxy counts 
(solid line) compared with a selection of literature data. In the 
case of the /iT-band counts we have followed the approach of 
ICristobal-Hornillos et aL t,2003i) and select only reliable counts 
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J vega 


0226 - 04 


2217 + 00 


1003 + 01 


1400 + 05 




N 


cr 




raw 


raw 


raw 


raw 


0.5 mag^'deg"^ 


0.5 mag^'deg"^ 




17.25 


28 


23 


26 


5 


758 


311 


63 


17.75 


32 


34 


30 


5 


934 


399 


36 


18.25 


59 


49 


42 


4 


1424 


770 


116 


18.75 


69 


56 


51 


13 


1748 


1087 


49 


19.25 


115 


80 


68 


21 


2627 


1889 


143 


19.75 


183 


129 


118 


33 


4283 


3510 


133 


20.25 


269 


169 


160 


39 


5892 


5056 


201 


20.75 


347 


254 


253 


45 


8316 


7462 


293 


21.25 


581 


350 


345 


65 


12404 


11534 


484 


21.75 


849 


470 


557 


99 


18268 


17365 


960 


22.25 


1183 


673 


969 


141 


27435 


26516 


2500 



Table 4. Columns 2 to 5 show, for each of our fields, the raw differential number counts per half magnitude bin. Column 6 show 
our raw total number densities (number O.Smag"' deg"^). Column 7 shows our total, differential galaxy densities, corrected for 
star contamination and Column 8 is their Icr error bar, obtained from the variance of the counts among the different fields. 



A'vcga 


0226 - 04 


2217 + 00 


1003 + 01 


1400 + 05 


N 

^ ' raw 


N 


cr 




raw 


raw 


raw 


raw 


0.5 mag^'deg"^ 


0.5 mag^'deg"^ 




15.75 


28 


17 


23 


4 


615 


274 


62 


16.25 


34 


29 


37 


5 


897 


467 


113 


16.75 


57 


51 


43 


14 


1409 


907 


136 


17.25 


84 


46 


62 


24 


1845 


1287 


89 


17.75 


128 


84 


100 


35 


2963 


2343 


42 


18.25 


228 


129 


136 


60 


4722 


4062 


284 


18.75 


287 


191 


217 


70 


6532 


5823 


284 


19.25 


406 


250 


313 


84 


8992 


8275 


469 


19.75 


545 


333 


394 


128 


11955 


11215 


453 


20.25 


836 


424 


543 


178 


16916 


16189 


780 


20.75 


1077 


516 


758 


261 


22304 


21596 


713 



Table 5. As in Tabled] but for iT-band counts. 



data from the literature, considering only data with negligible 
incompleteness correction and with star-galaxy separation ap- 
plied. Given the relatively small number of published 7-band 
counts we decided to plot most of the available data. It should 
be noted that we have been conservative in the selection of the 
magnitudes interval plotted in our counts, restricting ourselves 
to bins with relatively large numbers of galaxies, negligible in- 
completeness and small contamination corrections. The agree- 
ment with literature data is very good. For 7-band data the 
estimated slope of the galaxy counts in the range 17.25 < J < 
22.25, using a weighted least-squar es fit, is y.i ~ 0-39 ± 0.06, 
consistent with the findings of e.g. lSaracco et al.' ( 2001). The 
/T-band galaxy counts show an evident change of slope around 
K ~ 18.0. In the range 18 < K < 21.25 the slope of the 
galaxy counts jk ~ 0.29 + 0.08, with no significant hints 
of steeper slope down to the faintest magnitude levels. In the 
brighter magnitude range, 15.75 < TiT < 18, the slope is steeper: 
■Y k ~ 0.47 + 0.23 . Both results are consistent with the findings 
of ICiardner eTaP ((l996) and Cristobal-Hornillos et al. (2003 ). 
who also find a similar break, although at a slightly brighter 
magnitudes (K ~ 17.5). 



4.3. Star and galaxy colours 

In this section we will further evaluate the quality of our ab- 
solute and relative photometric calibration by investigating the 
colors of stars/galaxies in our fields. As the A'-band data have 
a slightly better seeing than /-band data (see Tables|2land|3}, 
to perform this analysis we used the sample of stars defined 
using Figure |6j 282 in total. We used SExtractor in dual- 
image mode to measure colours using matched apertures, us- 
ing /T-band images for detection and mag_AUTO, based on the 
/T-band flux distribution, to measure magnitudes in K and J 
bands. For the mosaicked images of field 0226 - 04, as usual, 
the chi-squared BVRIK was used as reference image for pho- 
tometric measurements, while stars were selected based on 
their /T-band image parameters. We checked that the differ- 
ence between colors obtained using MAG_AUTO are consistent 
with those one would obtain using a classical aperture magni- 
tude. Using only the good data quality area common to / and 
K-hand data, our final area totals 367 arcmin^. For the /-band 
data, for each field, we estimated as upper-limit for reliable 
magnitude measurements the one corresponding to a 3cr detec- 
tion within the circular aperture adopted by MAGJVUTO for faint 
and unresolved objects. Such upper-hmit values do not differ 
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Fig. 11. Differential 7-band number counts in our four fields. 
Error bars are Poissonian and no correction for stellar contam- 
ination has been applied to the points plotted. The dotted line 
shows the final total raw number densities, obtained by simply 
adding up raw densities from the four fields, while the heavy 
line shows the total galaxy densities obtained after correcting 
for star contamination. On the right of each points, slightly off- 
set for sake of clarity, are shown +1 - I cr error bars for total 
galaxy densities in each bin. See text for more details. 



significantly from those shown in Table|5]as 50% completeness 
values. For each field, 7-band magnitudes which were fainter 
than these values were replaced by the appropriate upper limits. 
Figure [21 shows the (7 - K) versus K colour-magnitude dia- 
grams for our data. The objects shown as star symbols indicate 
objects classified as point-like using our star/galaxy classifier 
(see section l4'.4.1> . 

It is reassuring to see that the majority of objects we clas- 
sified as stars are well separated from the global color distri- 
bution and are almost always below the dashed line (J - K) — 
0.96, corresponding to the typical color of a main sequence M6 
star. Furthermore the stellar locus is in same position for all 
four fields, indicating that our absolute calibration is accurate 
to within ~ 0.05 magnitudes. In Figure^] the dotted line at 
{J - K) = 2.3 corresponds to colour of an z ~ 2 evolved galaxy 
with a prominent 4000A break as a present day elliptical, or a 
z ~ 0.6 heavily reddened starburst galaxy. 

Figure^] shows galaxy (J - K) color distributions in dif- 
ferent bins of /iT-band magnitudes. It is evident how the popu- 
lation objects red in (7 - K) becomes progressively important 
at fainter magnitudes. A detailed analysis of the number counts 
and clustering properties of such red population will be pre- 
sented in a subsequent paper. 



Fig. 12. As in Figure^ but for /T-band. 
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Fig. 13. Our 7-band galaxy number counts, in units of num- 
ber 0.5mag"' deg"^ compared to a literature compilation, in- 
cluding counts from Tepli tz et al. ( 1999). S aracco et al. (^OOll) . 
iMaihara et allfcoOll) and lMartinilfcoOlT) . 



4.4. Clustering analysis for K- selected data 

In this Section we investigate the clustering properties of point 
like and extended sources in our /T-band catalogues. 

We use the projected two-point angular correlation func- 
tion, cl>{9), which measures the excess of pairs separated by 
an angle 6,9 + 66 with respect to a random distribution. This 
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Fig. 14. The ^T-band galaxy number counts, in units of 
number 0.5mag"' deg"^, obtained in this paper are com- 
pared to a cor npilation of those from the literature, includ - 
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Fig. 15. Colour-magnitude diagrams for the /T-selected sam- 
ples in each of our fields. The star symbols indicate objects 
classified as point-like using our star/galaxy classifier. See text 
for more details. 



17.5<K<18.0 
<J-K>=1.43±0.34 




Fig. 16. (J-K) colour distribution for the total /T-selected sam- 
ple in different slices of /T-band magnitude. Thjere is a trend 
towards redder colors at fainter magnitudes, probably reflecting 
an increasing fraction of high redshift, red, galaxies. 

statistic is useful for our purposes because it is particularly 
sensitive to any residual variations of the magnitude zero-point 
across our stacked imag es. We measure (jj{ff) using the standard 
iLandv & Szalavl(ll993h estimator, i.e.. 



DD - 2DR + RR 
RR 



(1) 



with the DD, DR and RR terms referring to the num- 
ber of data-data, data-random and random-random pairs be- 
tween and 6 + 56. We use logarithmically spaced bins, with 
Mogifi) - 0.2, and the angles are expressed in degrees, unless 
stated otherwise. DR and RR are obtained by populating the 
two-dimensional coordinate space corresponding to the differ- 
ent fields number of random points equal to the number of data 
points, a process repeated 1000 times to obtain stable mean val- 
ues of these two quantities. 

4.4.1. Clustering of point-like sources 

We first measure the angular correlation function ui ff) of the 
stellar sources. As stars are unclustered, we expect that, if our 
magnitude zero-points and detection thresholds are uniform 
over our field, then (jj(9) should be zero at all angular scales. 

The results for /T-band data are displayed in Figure [TTl 
where the correlation function is plotted for the total sample 
of stars obtained from all fields according to the procedure de- 
scribed in section : 282 stars in total for the /T-band images. 
At all scales displayed the measured correlation values are con- 
sistent with zero. Error bars are obtained through bootstrap re- 
sampling of the star sample (and are roughly twice poissonian 
error bars). A similar result is obtained for the measurement of 
clustering of stars obtained in 7-band images. 
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Fig. 17. Plot of the correlation function ll>{6) for the total sam- 
ple of stars in the /T-band stacks as a function of the logarithm 
of the angular pair separation in degrees. At all scales displayed 
the measured correlation values are consistent with zero. 

4.4.2. Clustering of extended sources 

The procedure followed to measure is similar to the one 
described above for the star sample. In the case of galaxies a 
positive amplitude of oj(ff) is expected, and we have to take into 
account the so called "integral constraint" bias. If the real ol>{0) 
is assumed to be of the form A^^ff', our estimator Q will be 
offset negatively from the true (jj{0), according to the formula: 



ojM = A^(6-' - C) 



(2) 



This bias increases as the area of observation decreases, and 
it is caused by the need to use the observed s ample itself to 
estimate its mean density, see eg'Pee ble'3 ( Il98d) . The negative 
offset AC can be estimated by doubly integrating the assumed 
true (jj(ff) over the field area Q: 



oj(0)d£lidQ.2 



(3) 



This integral can be solved numerically using randomly 
distributed points for each field: 



C 



(4) 



Assuming 1" as the pairs minimal scale at which two galax- 
ies can be distinguished as separated objects, and 6 = 0.8, 
we obtain the following values: C0226-04 = 7.58A, C22 17+00 = 
9.3 lA, C1003+01 = 8.60A, C1400+05 = 12.64A. 

We also measured the angular correlation function for our 4 
J- selected galaxy samples and verified that they do not exhibit 
any significant deviation from a power-law within the angular 
separation range associated to our sample areas. 



We estimated the amplitude A^ for a series of K limited 
galaxy samples by least square fitting A(0"'' ** - C) to the ob- 
served cl)(6), weighting each point using bootstrap error bars. 
Figure shows the results obtained for each of our fields on 
galaxy sub-samples of different /T-band limiting magnitudes. 
No correction for stellar contamination is applied (only the ob- 
jects classified as stars, using the method described in section 
14.4.11 were excluded from the analysis) and the error bars on 
the amplitude are Icr error bars obtained from rms of the fit. 
The different fields are in good agreement within the error bars. 
The 0226 - 04 field is the one with the larger area and the 
lower stellar contamination, and this probably explains why its 
amplitude values are systematically higher than those of the 
other fields. It should be remembered that in the presence of 
a randomly distributed spurious population among the sample 
of objects analized, like faint stars among our galaxy sample, 
the resulting measured correlation amplitudes are reduced by 
a factor (1 - /)^, where / is the fraction of the randomly dis- 
tributed component. Therefore a 10% (5%) contamination rate 
by stars implies a shifting of the values of the amplitude plot- 
ted by ~ 0.1 (0.06) downwards on the y-axis. Another point 
to consider is the expected cosmic variance among different 
fields. A simple numeric estimate of the expected variation of 
the measured amplit ude of the co rrelation function on the sky 
is o-A = aH^C^'^, see baddi etaP fcoOl ). For our fields the ex- 
pected scatter due to cosmic variance is roughly around 10% 
and comparable to the scatter observed among the fields. 

In Figure^]the filled circles show our final estimate for the 
amplitude of the correlation function. For each Umiting magni- 
tude these values are obtained by a weighted mean of the am- 
plitude of our four fields, and their error bars by computing the 
weighted variance among amplitudes (weighting proportional 
to the error bars on each field value). 

Figure shows the comparison of our results with liter- 
ature data. The con tinous, dotted and dashed lines show the 
models of PLE from lRoche et alJ il998'), with scaling from lo- 
cal galaxy clustering. A more detailed analysis, involving the 
use of spectroscopic (from the VVDS redshift survey) and pho- 
tometric redshifts for the galaxies of our sample is to be pre- 
sented in a forthcoming paper. 

Table|6llists our amplitude measurements for each field and 
for the total sample, in units of lO"'* at 1 deg, together with their 
Icr error bars, computed as explained above. For each field and 
limiting magnitude, the total number of objects N used in the 
analysis is also listed. 

5. Summary and conclusions 

In this paper we have presented a new near-infrared survey cov- 
ering four sub-areas located in each of the four fields of the 
VIMOS-VLT deep survey. We have described in detail our data 
reduction process starting from pre-reductions, astrometric and 
photometric calibrations, image resampling and stacking, and 
finally extraction of catalogues. At each stage, we have tried 
to quantify all sources of systematic and random errors in our 
survey. From extensive simulations, we have shown that our 
catalogues are reliable in all fields to at least K ~ 20.75 and 
J ~ 21.50: we define this limit as the magnitude where object 
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Table 6. Observed amplitudes A, in units of 10" at 1 deg, together with their Icr error bars, for each of our fields and for 
the total sample. For each field the number of objects used in the analysis down to the K magnitude faint limits shown in the first 
column is also listed. The number of objects in the first bin of field 1400 + 05 was too small, and therefore this bin was discarded 
from the analysis. 
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Fig. 18. Results obtained for the amplitude at 1 deg of each 
of our fields on galaxy sub-samples of different /T-band limit- 
ing magnitudes. No correction for stellar contamination is ap- 
plied and the error bars on the amplitude are Icr error bars ob- 
tained from rms of the fit. The different fields show a good 
agreement, within the error bars. On the right of each point, 
slightly offset for sake of clarity, filled circles show our final 
estimate, from the total sample, of the amplitude of the angular 
correlation function. 



contamination is less than 10% and completeness greater than 
90%. 

Based on repeated measurements of standard stars, we esti- 
mated that the error on our absolute photometric calibration, 
field-to-field, is ~ 0.05 magnitudes rm.s. This is consistent 
with the ~ 5% field-to-field variation of J- and /T-selected 
galaxy number counts and the measured field-to-field variation 
of the (J-K) colour of the stellar locus. We separated stars from 
galaxies using the parameter ri/2, which measures the radius 
for each object which encloses 50% of the total flux. Stellar 
counts for our four fields are consistent with the iRobin et al.i 



Fig. 19. We compare our total sample estimate of the ampli- 
tude, A^j at 1 deg, of the angular correlation function with a 
compilation of results from the literature. We include mesure- 
ments from: [Pad di et al. (2000), Kummel & Waaner (2001), 
McC racken et alJ ( l2000bl) . lDaddi et alJ ( 120031) and lRoche et all 
(I200.C The continous, dotted and dashed lines show, for refer- 
ence, the models of PLE from Roche et al. (1998) . with scaling 
from local galaxy clustering. 

(l2003h model of Milky Way, and our mean galaxy counts and 
counts slope over our four fields is also in excellent agreement 
with literature compilations. We observe a change in slope in 
the /T-band galaxy number counts at ~ 17.5 magnitudes. 

We investigated the colour-magnitude distribution of stars 
and galaxies identified in our catalogues. All objects lying in 
the K \s (J - K) stellar locus were successfully identified by 
our classifier For the galaxy population in the range 17.5 < 
K < 20.5 we measure a median (/ - K) colour of 1.4 + 0.3, 
consistent with published values. This value is remains approx- 
imately constant to progressively fainter magnitudes, until the 
faintest reliable limits of our sample (20.0 < K < 20.5). Our 
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fainter magnitude slices show some evidence of a red tail of 
objects (J - K) ~ 2.5 which becomes progressively larger at 
fainter magnitudes. 

Finally, we measure the angular clustering of stars and 
galaxies for our four fields. Our stellar correlation function is 
consistent with zero for all four fields on all angular scales. 
The amplitude of our galaxy correlation function shows the 
expected scaUng behaviour for increasingly fainter magnitude 
slices, and is consistent previously-presented measurements. 

These catalogues will be an excellent tool to investigate the 
properties of distant galaxies selected in the near-infrared, and 
such investigations which will be the subject of several forth- 
coming articles. 
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